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Abstract 

It was recently found that the electric local-field effect (LFE) can lead to strong coupling of atomic Bose-Einstein condensates 
(BECs) to off-resonant optical fields. We demonstrate that the magnetic LFE gives rise to a previously unexplored mechanism 
for coupling a (pseudo)spinor BEG or fermion gas to microwaves (MWs). We present a theory for the magnetic LFE, and 
hnd that it gives rise to a short-range attractive interaction between two components of the (pseudo) spinor, and a long-range 
interaction between them. The latter interaction, resulting from deformation of the magnetic held, is locally repulsive but 
globally attractive, in sharp contrast with its counterpart for the optical LFE, produced by phase modulation of the electric 
held. Our analytical results, conhrmed by the numerical computations, show that the long-range interaction gives rise to 
modulational instability of the spatially uniform state, and creates stable ground states in the form of hybrid matter-wave- 
microwave solitons (which seem like one-dimensional magnetic monopoles), with a size much smaller than the MW wavelength, 
even in the presence of arbitrarily strong contact inter-component repulsion. The setting is somewhat similar to exciton- 
polaritonic condensates in semiconductor microcavities. The release of matter waves from the soliton may be used for the 
realization of an atom laser. The analysis also applies to molecular BECs with rotational states coupled by the electric MW 
held. 

PACS numbers: 03.75.Lm, 05.45.Yv, 42.65.Tg 


Ultracold atomic gases are used in a various areas, 
including quantum metrology and interferometry i-i, 
and the emulation of nonequilibrium ^antum dynamics 
0 and condensed-matter physics Q-S- They have also 
drawn much interest as tunable media for quantum op¬ 
tics. In this vein, co-manipulation of quantum light and 
matter waves has been studied in cavities loaded with 
atomic Bose-Einstein condensates (BECs) @, 113 . Ra¬ 
man superradiance in ultracold gases trapped in a cavity 
was used to generate stationary lasing with a bandwidth 

< 1 MHz, and with the average cavity photon number 

< 1 0. A mirrorless parametric resonance has been 
demonstrated for atomic BEC loaded into an optical lat¬ 
tice (OL) [ 13 . Optomechanics-induced large-scale struc¬ 
turing of ultracold atomic gases has been reported in Ref. 
[ 13 . The resonant interaction of laser fields with BEC 
was also proposed for generating “photonic bubbles” em¬ 
ulating cosmology settings [l3| . 

An important feature of the interaction of light with 
ultracold gases is the local-field effect (LEE), i.e., a feed¬ 
back of the BEC on the light propagation. Strong LEE 
can be induced in cold-atom experiments, as recently 
demonstrated with the help of OLs [iM3- Usually, OLs 
are sturdy structures, maintaining perfect interference 
fringes. However, asymmetric matter-wave diffraction 
on an OL formed by counterpropagating optical fields 
with unequal intensities [l3 could be explained only by 
taking into regard deformation of the OL by the LEE 
[ 13 . Conventional rigid OLs and their deformable coun¬ 


terparts may be categorized as “stiff’ and “soft” ones. 
Polaritonic solitons, produced by hybridization of cou¬ 
pled atomic and optical waves, have been predicted in 
soft OLs [13 , being promising for the matter-wave inter¬ 
ferometry due to their high density [13. These results 
demonstrate the potential of the soft OLs in studies of 
systems combining quantum matter and photons, akin 
to exciton-polaritons in microcavities [^ . Furthermore, 
BECs built of up to 10® atoms are now available [2lj. For 
such massive BECs, the refraction-index change through 
the perturbation of the atomic density may be significant, 
even for the laser-frequency detuning from the resonance 
^ 1 CHz, allowing the LFE to generate hybrid matter- 
wave-photonic states El- 

The use of spinor gases opens ways for the emulation of 
the spin-orbit coupling [3 and quantum magnetism [23 , 
as well as for the realization of quantum matter-wave [23l| 
and microwave (MW) [13 optics. In this context, the 
MW magnetic field is used for manipulating spin states. 
Coupling different hyperfine atomic states by MWs was 
studied in other contexts too, includin g d ressed states 
[13 , domain walls [13 , and instabilities (27j | . 

However, manifestations of the magnetic LFE (MLFE) 
in quantum gases were not studied yet, unlike its electric 
counterpart. In this work, we develop the theory of the 
MLFE for a MW field coupled to the pseudospinor BEC. 
The MW wavelength (> several mm) exceeds the typ¬ 
ical size of the BEC by orders of magnitude. In this 
situation, the BEC was considered before as a thin slice 
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that affects the phase of the MW field. We find that the 
MLFE causes subwavelength deformations of the MW 
amplitude profile too, inducing a long-range interaction 
between components of the pseudospinor BEC. Unlike 
the electric LEE , where nonlocal interaction is in¬ 
duced by phase perturbations, the long-range interaction 
generated by the MLFE is locally repulsive but glob¬ 
ally attractive. The same effect leads to local attrac¬ 
tion between the components of the BEC, which may 
compete with collisional repulsion between them. We 
demonstrate that these interactions create self-trapped 
ground states (GSs) in the form of hybrid matter-wave- 
microwave solitons, whose field component seems like 
that of a magnetic monopole. Actually, the solitons 
realize a dissipation- and pump-free counterpart of hy¬ 
bridized exciton-polariton complexes in dissipative mi¬ 
crocavities, pumped by external laser fields. Opposite to 
our case, the size of those complexes is much larger than 
the polaritonic wavelength, while the effective mass of 
the excitons is usually considered infinite [l^. Note that 
direct dissipation-free emulation of the exciton-polariton 
setting is possible too in a dual-core optical system 
We consider the magnetic coupling of the MW radi¬ 
ation with frequency wl to two hyperfine atomic states 
II) and It) which compose the BEC pseudospinor @,[11, 
with free Hamiltonian Hq = fp' /{2m) — (M/2) 0 - 3 , where 
M is the energy difference between the two states, (T 3 
is the Pauli matrix, and p the atomic momentum. The 
magnetic interaction is governed by term 

Hint = - ■ (Be““^‘ -h . (1) 

Here, are matrix elements of the magnetic 

momentum, and the magnetic induction is B = 
/TqH -|- M with M = with 1 /^ and 

building the pseudospinor wave function, l^h) = 

exp 6 xp (—MLt/2)) . In the rotating- 

wave approximation, ly/) = satisfies the system 

of coupled Gross-Pitaevskii equations (GPEs): 


Helmholtz equation to the normalized form. 
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' dr 
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(3) 

din + tPn = ( 4 ) 

with X = X/Aq, 

n=H/hc and /3 = 1114 ,-^ • ni-^ito/{fiSXo),ii = toA/2, 
7 = ]>leoUjl^mifXQ/{hcS), k = Xqlol/c measured 
in natural units of time and coordinate, to = 
h/{fj,Qmi^hc), Xq = ^Jhto/m, where he is a magnetic- 
field strength and N the total atom number. The rescaled 
wave function is subject to normalization = 


f-x 1 ^ 4 .(^)P] = 1 - If collisions between 

atoms in different spin states are taken into account, 
which may be controlled by the Feshbach resonance [11 , 
P in Eq. m combines contributions from the MLFE 
and direct interactions. On proper rescaling, the same 
system of Eqs. ®,0 applies to a degenerate gas of 
fermions [SlL l32l| with spin 1/2, in which tjj^ and ijj^ rep¬ 
resent two spin components, coupled by magnetic field 
n, and asymmetry rj is imposed by dc magnetic field. 

Equati on pi ) can be solved using the respective Green’s 
function [33|, n{x,t) = Aexp{iKx) + C exp{—iKx) — 


( 7 / 2 k) sin (k \x — a;'|) (/J (s', t) {x', t) dx', where 
constants A and C represent the solution of the corre¬ 
sponding homogeneous equation, i.e., they are ampli¬ 
tudes of two incident counterpropagating microwaves. 
Since the MW wavelength is far larger than the con¬ 
densate size (k ^ 10 “^), one may set exp(zKx) ~ 
1 and s\'a.{K\x — x'\) /K ~ \x — x'\ in the domain oc¬ 
cupied by the condensate, to simplify the solution: 


n{xn) = m-{i/2)C^ 


x — X 


(j)’/ {x', t) (j)^ {x', t) dx', 


where "Hq = A-l-C is made real by means of a phase shift. 
The form of the magnetic field at |a:| —>■ 00 resembles that 
of a one-dimensional artificial magnetic monopole [sl : 
"Hasympt (a:) = -jN^lx/ (l>l (x) (j)y (a;) dx. 

Substituting the solution for % in Eq. , we arrive 
at the final form of the GPEs: 
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1 d'^cj), I ,2 

(j)^ {x',t) ipl {x',t) dx' , (5) 


with the MW detuning from the atomic transition A = 
ujl — 5. Neglecting the time derivatives of H and M for 
the low-frequency MW, the wave equation for H reduces 
to the Helmholtz form, V^H -|- w^/c^H = —eow|M. 

We consider a cigar-shaped condensate with effective 
cross-section area S, subject, as usual, to tight transverse 
confinement, and is irradiated by two linearly-polarized 
counterpropagating microwaves along the cigar’s axis x. 
Eliminating the transverse variation of the fields under 
these conditions [ 2 ^, we reduce the coupled GPEs and 


dd).. 1 I ,2 

/ + 00 

\x — x'\(j)*^{x',t) 4 >.^{x',t) dx'. ( 6 ) 

-00 

Thus, the MLFE gives rise to two nonlinear terms: the 
one ^ P accounts for short-range interaction, while the 
integral term represents the long-range interaction, which 
is locally repulsive, but globally attractive because the re¬ 
pulsion kernel, |a; —a:'|, growing at |a;| —>■ 00 , suggests 
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a possibility of self-trapping. The mechanism of creat¬ 
ing bright solitons by the spatially growing strength of 
local self-repulsion was proposed in Ref. [35|, and then 
extended for nonlocal dipolar-BEC [s^ settings; how¬ 
ever, that mechanism was imposed by appropriately engi¬ 
neered spatial modulation of the nonlinearity, while here 
we consider the self-trapping in free space. 

The symmetric version of Eqs. m and dH), with 
ry = 0 , may be combined into separate equations for 
4’± = 4>i ^ with trapping (for -|-) and expulsive (for 
—) potentials, respectively. Therefore, this system has 
only symmetric solutions {(j)_ =0). At \x\ —>■ oo, Eqs. 

and ([ 6 ]) take the linear asymptotic form, = 

— + ( 7 A^j:/ 2 ) hence solutions for the 

asymmetric system (ry 0 ) have symmetric asymptotic 
tails too, (j)^ = (j)^ ^ exp (—( 2 / 3 )y' 7 dV^|a;p/^). 

The GS of system ([5]), dH) with 7 y = 0 and chemical 
potential fi is sought for as (j)^^ = where real 

ip satisfies equation 


flip = 


2 dx"^ 




\x — x'\ {x') dx' 




(7) 

with fi = fi — Hq. Thus, Hq only shifts the chem¬ 
ical potential in the zero-detuning system. For /3 = 
0, it follows from Eq. © that fi and the magnetic 
field obey scaling relations: {//(y),"H (a:; 7 ) — "Ho} = 

( 7 / 70 )^^^ {d (7 = 7o). 7^ ( 2 ^; 7 = 7o) - 77o}, where 70 is 
a fixed constant. Thus, for ?y = 0 and /3 = 0, all the GSs 
may be represented by a single one, plotted in Fig. [Ija), 
which was found by means of the imaginary-time method 
. This is a hybrid soliton, built of the self-trapped mat¬ 
ter wave coupled to the subwavelength deformation of 
the magnetic field. 

In the presence of the local self-repulsion (/3 < 0), the 
GS can be found with the help of the Thomas-Fermi ap¬ 
proximation (TEA), which neglects the second derivative 
in Eq. dZ]): 




(1/4)\/7/|/3|cos^, at |^| < tt/2, 
0, at lei > 7 r/ 2 , 


"Htfa (a:) + fi 


(1/4) cos C, at 1^1 < 7r/2, 

(l?l - 7 r/ 2 ) /4, at |^| > 7 r/ 2 , 

(9) 


where ^ = \J^/\j3\x. An example, displayed in Fig. [Ijb), 
shows very good agreement of the TEA with the numeri¬ 
cal solution. Thus, the globally attractive long-range in¬ 
teraction induced by the MLFE creates the self-trapped 
GS, overcoming the arbitrarily strong self-repulsive con¬ 
tact interaction. 

The time-of-flight spectrum produced by releasing the 
condensate can be used in the experiment to detect 
the solitons predicted here, characterized by their dis¬ 
tribution over the longitudinal momentum, (j) {kx) = 
e~'‘^^^p{x)dx, as shown in Fig. [T] for /3 = 0 and 



FIG. 1: (Color online) (a) The GS wave function, p{x), along 
with its Fourier transform, (j}{kx), and magnetic field, H (x), 
for d ~ 0, with p = 4.26. (b) Comparison of the GS, as 

predicted by the TFA [Eqs. dS)! , dll) and (IT0I1 ] (dashed curves) 
and found numerically (solid lines) for /3 = —5 with p — 0.079. 
In both plots, 7 = 0, Ho = 0, and 7 = 10~®. 


—5. In particular, the TFA yields 


4’tfa (kx) — 


•(1/31/7) 


1/4 


ST (5/4 + ^/Wl^kx) r (5/4 - ^/WJ^kx) 

( 10 ) 

where T is the Gamma-function, see the right bottom 
panel in Fig. [TJ Note that expression m vanishes at 
kx = ± 2 y^ 7 /|, 5 | (5/4 -I- n ), n = 0,1,2,... . The strong 
compression of the soliton in the momentum space, evi¬ 
dent in Fig. 
wave lasers 


TJR), may be used for the design of matter- 
37l |. as the released beam will feature high 
velocity coherence. 

The existence of bright solitons in free space is related 
to the modulational instability (MI) of flat states [ 3 ^ . 
The flat solution to Eq. © is (/ = $0 exp {—ip^t ), 
with divergence of p^ regularized by temporarily re¬ 
placing \x — x'\ in Eq. dZl by |x — a;'| exp (—e [x — a:'|) 
with small e > 0. The MI analysis for small pertur¬ 
bations with wavenumber k and MI gain A proceeds, 
as usual, by substituting = $ (x, t) exp (fy (x, t)), 

with {$;x} = {5 >o;-MoO + {^l°^X®}exp(ifca:-h At), 

where are perturbation amplitudes. The 

subsequent linearization and then setting e = 0 yields 
A^ = — — /3 $qA:^ — 7 $q) . Due to the nonlocality, 

A^ does not vanish at —>■ 0, in contrast with local mod¬ 
els [ 3 ^. The MI is always present, as A^ remains positive 

at < 2 1$ 


oil V P '^0 + 7 + /3 l^olj ■ Thus, arbitrarily 

strong local self-repulsion, with /3 < 0 , does not suppress 
the MI. 

In the system with detuning, i.e., 7 7 ^ 0 in Eqs. dS]) 


3 


























FIG. 2: (Color online) Numerically found profiles of the GS 
wave functions and magnetic field in the system with rj = 1, 
/3 = 0, and 7 = 10“®, for Ho = 0.5 (a), 1.5 (b), 6.0 (c). The 
respective chemical potentials are fj, = —1.12 (a), —1.80 (b), 
—6.079 (c). For the strongly asymmetric soliton in (a), the 
analytical approximation predicts the amplitude ratio of the 
two components 0.250, while its numerically found counter¬ 
part is 0.235. 



FIG. 3: (Color online) The relative share of the total norm of 
each component in the system with 7 = 10“®, P = 0, rj — 1, 
vs. the background magnetic field, Ho. The inset shows the 
dependence for ^4. at small values of Ho (stars) vis-a-vis the 
analytical prediction. 


and ©, the background magnetic field Hq is an essential 
parameter. Figure[5]plots the GS wave functions, {x), 
and the corresponding magnetic field, H (x), for different 
values of Ho. The GS exhibits asymmetry between the 
lower- and higher-energy components at Ho ^ ry, while 
large Ho suppresses the asymmetry, as seen in Fig. [31 
which displays the scaled norms of the two components, 
versus Hq. 

For ?7 ^ Ho, strongly asymmetric GSs can be found 
using the stationary version of Eqs. m and (0) with 
chemical potential p. = —ry -|- Ayr, |A/r| <C rj. Then, Eq. 
m eliminates the weak component in favor of the strong 
one: 75ft! {2r])~^ and Ay ft Hq/ (4ry^), which 

agrees well with numerical results, as shown by the inset 


in Fig. [3] The substitution of this into Eq. ([6]) yields 

1 /3H2 


Ayt 


27? 


2 dx^ 

'+00 


4?y- 


0 


0/2 y+oo 

+i^n(^)y_ \x-x'\p^x')dx', (11) 


which is actually tantamount to Eq. ©• The respective 
small deformation of the magnetic field is H (x) = Ho — 
(7Ho/4?y) (j)\{x')\x-x'\ dx'. 

We have confirmed the stability of all the GS states 
by direct simulations of Eqs. (O and (jS]) with randomly 
perturbed initial conditions. For the symmetric system 
with 7y = /3 = 0, the above-mentioned scaling implies that 
the stability of a single GS guarantees the stability of all 
GSs, while for the detuned system the stability had to 
be checked by varying Ho at fixed 7 and ?y. Further, the 
stability for ?? = /3 = 0 is predicted by the anti- Vakhitov- 
Kolokolov (anti-VK) criterion, which states that the nec¬ 
essary stability condition for bright solitons supported by 
the repulsive nonlinearity is dp/dN > 0 (the VK cri¬ 
terion proper, which pertains to attractive nonlinearity, 
is dp/dN < 0 [3^,|4l|). Although N = 1 was fixed above, 
the criterion can be applied by means of rescaling which 
fixes 7 and liberates N. As a result, the scaling relation 
p 7^/3 is replaced by yt hence the criterion 

holds. 

Summarizing, we have explored the MLFE (magnetic 
local-field effect) in the BEG built of two atomic states 
coupled by the MW (microwave) field. We have deduced 
the system of evolution equations for the matter-wave 
components and MW magnetic field, which demonstrate 
that the subwavelength distortion of the magnetic field 
by perturbations of the local atom density induces short- 
and long-range interactions between the BEG compo¬ 
nents. The same equations apply to the spinor wave func¬ 
tion of a fermionic gas coupled to the MW magnetic field. 
The model produces the self-trapped GS (ground state) 
in the form of the hybridized BEG-MW subwavelength 
solitons, which may be considered as counterparts of hy¬ 
brid exciton-polariton solitons in the dissipation-free sys¬ 
tem. Basic characteristics of the solitons were obtained 
analytically. The release of the solitons from the cigar¬ 
shaped trap may be used as a source of coherent matter 
waves for an atom laser. The flat states in the present 
system are subject to the modulational instability, which 
is naturally related to the existence of the bright solitons. 

It is straightforward to extend the analysis for molec¬ 
ular BECs, with the transition between two rotational 
states driven by the electric MW field. Such ultracold 
molecular gases have a potential for quantum simula¬ 
tions of condensed-matter physics [i^. An interesting 
extension may be the analysis of the system with a three- 
component bosonic wave function corresponding to spin 
H = 1 [4^, in which a single MW field couples compo¬ 
nents with mp = ±1 to the one with mp = 0- Other 
relevant directions for the extension are search for ex¬ 
cited states in the system, in addition to the GS, and the 
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analysis of the two-dimensional setting. Furthermore, it 
is relevant to investigate a potential effect of the MLFE 
on quantum precision measurements. 
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